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A comprehensive analysis of optical second harmonic generation (SHG) on an incommensurate (IC) magnet- 
ically ordered state is presented using multiferroic MnW04 as model compound. Two fundamentally different 
SHG contributions coupling to the primary IC magnetic order or to secondary commensurate projections of 
the IC state, respectively, are distinguished. Whereas the latter can be described within the formalism of the 
122 commensurate magnetic point groups the former involves a breakdown of the conventional macroscopic 
symmetry analysis because of its sensitivity to the lower symmetry of the local environment in a crystal lat- 
tice. Our analysis thus foreshadows the fusion of the hitherto disjunct fields of nonlinear optics and IC order in 
condensed-matter systems. 
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I. INTRODUCTION: IMAGING INCOMMENSURATE 
ORDER 

A crystal structure is called incommensurate (IC), if a pe- 
riodic deviation of atomic parameters from the basic crys- 
tal structure cannot be related to the periodicity of the un- 
derlying crystal structure by a rational number^ The viola- 
tion of the translation symmetry of the basic crystal struc- 
ture by the IC long-range order manifests in an astonishing 
diversity of unusual physical effects so that systems with IC 
phases attracted a lot of attention since Dehlinger reported the 
first aperiodic crystal in 1927,- Nowadays, IC order is known 
to play an essential role in quasicrystals^, liquid crystals^i 
magnetic multilayers,— or polymers.— Other examples for the 
omnipresence of IC phases can be found among systems 
with strong electronic correlations such as high-temperature 
superconductors)^ colossal-magnetoresistance manganiteSj^ 
or multiferroic Si^ The pronounced reduction of symmetry by 
the IC order reflects the complexity of the underlying micro- 
scopic interactions. Macroscopically, this can lead to a variety 
of phases that become allowed in the IC statai^ and enrich the 
phase diagram of the host material — a key ingredient for the 
design of multifunctional materials. 

Despite the longstanding interdisciplinary interest in IC 
structures, their theoretical and experimental analysis is still 
a challenge. For instance, it was realized only recently that an 
IC magnetic spiral may induce a spontaneous electric polar- 
ization, and a controversial discussion about the microscopic 
mechanisms of this coupling is persisting jii"— For analyzing 
the order parameters and domain structures of an IC struc- 
ture, a group-theoretical approach based on Landau theory 
of phase transitions is usually chosen. A description of the 
IC state in a superspace with > 3 dimensions'^ or, alterna- 
tively, by considering the irreducible representations of the 
three-dimensional space group in the presence of a modulated 
structure^ are possible. Experimentally, neutron and X-ray 



diffraction are usually employed for investigating aperiodic 
spin, charge, or lattice modulations. Yet, although these es- 
tablished techniques constitute a versatile tool for probing IC 
structures, essential ingredients characterizing the IC state re- 
main inaccessible. In particular, this applies to the domains 
associated to the IC state. On the one hand, the spatial resolu- 
tion of a diffraction experiment is often limited because of the 
integrating nature of the experiment. On the other hand, most 
of the IC magnetic structures do not carry a macroscopic mag- 
netization so that it is difficult to couple to the order parameter 
and its spatial distribution. However, since domains and do- 
main walls determine many macroscopic physical properties 
and, thus, the technological feasibility of a compound, access 
to these domains in a convenient way is highly desirable, par- 
ticularly because it is expected that the IC nature of the or- 
dered state will lead to domains with novel properties. 

A technique with inherent spatial resolution and access 
even to "hidden" ordered structures is optical second har- 
monic generation (SHG)J^ For decades SHG has been ap- 
plied to image structural, magnetic, or electric domains and to 
reveal spatial correlation effects between them .'^i'^ More re- 
cently, SHG was even used to resolve the sub-picosecond dy- 
namics of a magnetic order parameter . '^i^" However, almost 
any form of long-range order probed by SHG thus far was 
commensurate. In the very small number of investigations on 
IC structures it was demonstrated that the emergence of the IC 
state is reflected by changes in the SHG yield, but a theoret- 
ical analysis of the detected signal was never attempted.—"— 
The same holds for recent investigations on IC magnetic ferro- 
electrics where SHG was used for imaging magnetic and fer- 
roelectric domains without performing a detailed discussion 
of the SHG contributions and the IC symmetries involved . 
Thus, a resilient framework categorizing the nonlinear inter- 
actions between the electromagnetic light fields and an aperi- 
odically ordered state of matter is yet to be developed. 

In this report we present a comprehensive analysis of opti- 
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cal SHG on an IC magnetically ordered state. Using multifer- 
roic MnW04 as model compound, two principally different 
SHG contributions coupling to the symmetry-breaking order 
parameters are identified. The first one corresponds to a com- 
mensurate projection of the IC magnetic state that is repre- 
sented by an improper order parameter and can be described 
within the 1 22 commensurate magnetic point groups. The sec- 
ond one reflects the entire IC structure of the compound and 
involves a breakdown of the conventional macroscopic sym- 
metry analysis since local symmetries become involved. With 
our analysis we foreshadow the fusion of the hitherto disjunct 
fields of nonlinear optics and IC order in condensed-matter 
systems. 
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FIG. 1: Overview of the magnetic phase a • . the symmetry- 
breaking order parameters^^iS, and the corresponding point groups 
for muUiferroic MnW04. The commensurate AFl phase at T < 
7.6 K is not discussed in this manuscript and therefore not shown. 



II. INCOMMENSURATE ORDER IN MULTIFERROIC 
MnW04 

Multiferroics are compounds uniting at least two primary 
forms of (anti-) ferroic order in the same phase of a mate- 
rial. The majority of research activities is devoted to the mag- 
netically ordered ferroelectric s since the coexistence of mag- 
netic and electric order is a potential source for magnetoelec- 
tric cross-correlations allowing one to manipulate magnetic 
order by electric fields (or vice versa). The most pronounced, 
so-called "gigantic", magnetoelectric coupUng effects are ob- 
served when IC magnetic order promotes the emergence of 
a spontaneous electric polarization that is rigidly coupled to 
the magnetic order parameter. Because of the novel and un- 
usual nature of this manifestation of IC order and the intense 
desire to advance our understanding of it, the choice of an IC 
magnetically induced ferroelectric as model system for unrav- 
elling the relation between IC order and SHG suggests itself. 

We select MnW04 because it combines a relatively sim- 
ple crystallograhic structure with a variety of IC phases. The 
magnetic phase diagram is well established and the gen- 
eral feasibility of an approach by SHG has already been 
demonstrated^^"— so that we can now focus on analyzing 
the correlation between SHG and the IC state. MnW04 is 
closely related to compounds such as TbMnOa, Ni2V308, or 
CuO^"— so that the results gained here will be of general use. 

The multiferroic phase of MnW04 is governed by two 
fundamentally different types of order parameters: (i) a set 
of two primary magnetic order parameters (tJaf,, ^Tafz), and 
(ii) one secondary electric order parameter {P^F). The se- 
quence of phase transitions leading to the two classes of or- 
der parameters in MnW04 is sketched in Fig. [T] Below 
Tn = 13.5 K the magnetic moments of Mn^+ align in a 
collinear way along the easy axis while their magnitudes 
are sinusoidally modulated with k^Fj = (—0.214,5,0.457)^. 
The resulting IC spin-density wave breaks the translation 
symmetry of the lattice but conserves the inversion symme- 
try of the crystal. This so-called AF3 phase can be de- 
scribed by a single two-dimensional magnetic order parame- 
ter T7af3 = (o'AFje'^'^-'jffAFje"'^'^')^^. Upon further cooling an 
additional transverse spin component orders at T2 = 12.7 K. 
In this AF3 — !> AF2 transition the spin-density wave turns 
into an elliptical spin spiral without changing the IC periodic- 



ity of the modulation so that k^F^ — k^F^- The AF3 ^ AF2 
transition is driven by a second magnetic order parameter 
riAfi ^ (dAFze'^*''-, CJAF2e"'^*''^)- The coexistence of 77af2 and 77af3 
breaks the inversion symmetry and induces a spontaneous po- 
larization Py^ along the y axis {y\\b) according to ^^'^^i^^ 

i'f - CJaf.(Taf3 , (1) 

thus constituting the multiferroic phase of MnW04. The cor- 
relation expressed by equation ([T]) implies that the IC nature 
of the primary (proper) order parameters t], with amplitudes 
(T; is projected onto the secondary (improper; here: pseudo- 
proper—) order parameter Py^. The projection corresponds 
to a commensurate net polarization that may be additionally 
modulated in an IC way around its net value^^ Note that 
the incommensurability of Py'' is not mandatory and, even 
if present, remains mostly unnoticed because experimental 
techniques probing the spontaneous polarization such as py- 
roelectric measurements measure the commensurate contribu- 
tion only. Macroscopically the polarized state can be therefore 
described by the point-group formalism used for commensu- 
rate structures, as in the lower part of Fig.[T] 

III. OPTICAL SECOND HARMONIC GENERATION 

An experimental method that is known to be particularly 
sensitive to the point symmetry of a crystal is optical SHGi^. 
It is described by the equation^ 

Si{2co)^eoX,jkEj{(o)E,,{co) . (2) 

An electromagnetic light field E at frequency (o is incident 
on a crystal, driving a charge oscillation S(2a)), which acts 
as source of a frequency-doubled light wave of the intensity 
^SHG °^ |S(2a))p. The source term S can be written as multi- 
pole expansion 

, d^P(2co) / dM(2co)\ fd^(VQ)\ 

(3) 

with P(2co), M(2a)), and VQ{2co) as electric -dipole (ED), 
magnetic-dipole (MD), and electric -quadrupole (EQ) contri- 
bution, respectively. The nonlinear susceptibility Xijk cou- 
ples incident light fields with polarizations j and A; to a SHG 



3 



Point 


focident MD/EQ-SHG 


ED-SHG 


group 


light 


contributions 


contributions 


lylmyV 


K 1 X 


,yMD.EQ yMD,EQ ,yMD,EQ 

Avvv ' ^yzz ' Azyc 
















yMD,EQ yMD,EQ yMD.EQ 




2,1' 




yMD,EQ yMD.EQ yMD,EQ 

Avvv ' '^yzz ' ^zyz 


yED yED yED 

A_vv_v" ^yzz'' ^zyz 


(AF2) 




yMD,EQ yMD.EQ yMD,EQ 


yED yED yED 

XvVV A-V.V.V' >C.VV.V 



TABLE I: Non-zero SHG tensor components Xiji^ allowed on the ba- 
sis of a symmetry analysis within the 122 magneto-crystalline point 
groupsi^ The SHG tensor is related to a Cartesian reference system 
(yX, y, z) as defined in [3^ Calculations were done for the magnetic 
AF3 phase and the multiferroic AF2 phase of MnW04 for light prop- 
agating along the x, y, and z axis. Note that SHG contributions with 
! 7^ i cannot be detected with light propagating along x, y, or z 
so that they are omitted. MD-SHG and EQ-SHG lead to the same 
selection rules and are therefore not distinguished. 



wave with polarization ;. According to Neumann's principle 
the symmetry of a compound determines the set of nonzero 
components As a consequence, ED-SHG vanishes in 

centrosymmetric media while MD-SHG and EQ-SHG remain 
allowed. 

Using Eqs. (|2]) and Q and applying symmetry-dependent 
selection rules^ the tensor components summarized in Ta- 
ble H] are derived as contributions to SHG in MnW04. Note 
that TableUis based on the use of point groups (as in Fig.[T]i- 
This is based on the assumption that translation operations 
and translation symmetries can be neglected in optical exper- 
iments since they are not resolved by the light. More quanti- 
tatively, the variation of the electromagnetic amplitude of the 
light field across the expansion of a unit cell (as limit of the 
non-primitive lattice translations) is considered to be small 
enough to be neglected. With respect to Table H] this means 
that the aperiodic nature of the order parameters with their 
violation of the translation symmetry is also not taken into ac- 
count. The corresponding SHG light will be termed commen- 
surate (C-) SHG in the following in order to emphasize that 
averaged, lattice-periodic properties are considered. Table U 
reveals that for light fields propagating along the y axis (A: || y) 
SHG is forbidden in all phases. For | [ a or 1 1 z MD-SHG and 
EQ-SHG can occur at any temperature whereas ED-SHG is 
restricted to the multiferroic AF2 phase. 

A detailed discussion of the technical aspects of SHG in 
ferroic systems including the experimental setup used here is 
found in Refs. [l^ and |27l In the present experiment, melt- 
grown MnW04 single crystals^ were processed into sets of 
polished platelets with a thickness of about 100 /im with each 
of the samples being oriented perpendicular to the one of 
the Cartesian axes x, y, zM- The samples were illuminated at 
normal incidence in a transmission setup by light pulses of 
2 — 5 mJ and 3 — 8 ns at a repetition rate of 10 — 40 Hz. 
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FIG. 2: Optical transitions and paramagnetic SHG contributions in 
MnW04. (a) Linear absorption spectrum of MnW04 obtained with 
A-polarized light propagating along the y axis of the crystal. At low 
temperature two intra-atomic and one charge-transfer transitions are 
identified and labelled accordingly, (b) SHG spectra obtained in the 
paramagnetic phase with light incident along the x axis. Only crys- 
tallographic SHG from x"^,^ is observed. 



IV. EXPERIMENTAL RESULTS AND DISCUSSION 

Based on Table U the experiments will proceed as follows: 
A. identification of electronic transitions in the optical range 
and of any SHG contributions not coupling to magnetic or 
electric long-range order; B. identification and analysis of C- 
SHG contributions; C. identification and analysis of IC-SHG 
contributions. 



A. Optical transitions and paramagnetic SHG contributions 

For identifying the electronic transitions of the Mn^+ ion 
and the crystallographic background contributions to SHG, 
Fig. |2] shows linear absorption spectra and polarization- 
dependent SHG spectra in the paramagnetic phase of 
MnW04. The absorption spectrum taken with x-polarized 
light incident along the y axis exhibits a steep increase beyond 
2.7 eV which corresponds to the lowest 0^^-Mn^+ charge 
transfer .^^1^^ In addition, absorption peaks with a width in the 
order of 0. 1 eV are observed and marked by black arrows in 
Fig. 121 a). These peaks at 2.18 eV and 2.65 eV are assigned 



to the intra-atomic 



'Tu and M 



T20 transitions 



between the Mn^+(3t/^) orbitalsi^ 

In agreement with the symmetry analysis ED-SHG contri- 
butions were not observed above 7n. The signal in Fig.|2b) 
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FIG. 3: Commensurate SHG in the multiferroic phase of MnW04 
(arbitrary units), (a, b) C-SHG spectra obtained with k\\x (Xy?z' Xy°y) 
and k\\z (Xrxx)- The exceptionally large signal yield at 1.95 eV due 
to multiferroic phase-matching, (c) Temperature dependence of the 
SHG intensity from x^,, X\°v' ™d X\°x- '^e signals are associated 
to ED-SHG contributions that are present in the multiferroic AF2 
phase only. 
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FIG. 4: Incommensurate SHG in the magnetically ordered phases of 
MnW04. (a) IC-SHG spectra obtained with yt|| v for which C-SHG 
contributions of any multipole order are not allowed. The temper- 
ature dependence in the inset confirms the presence of the IC-SHG 
signal in the multiferroic AF3 phase and the non-multiferroic AF2 
phase, (b, c) Polarization dependence of the IC-SHG signal at 8 K 
and 2.22 eV for light incident parallel to the y axis, x' and z' de- 
note the axes of the rotated coordinate system consistent with the 
polarization dependence of the IC-SHG signal. P and A denote the 
polarization of the incident light and of the detected SHG compo- 
nent, respectively. In (b) P and A are rotated simultaneously while in 
(c) P is rotated while A remains constant. 



shows the ^Aig — > ^Tig transition and is attributed to the J™'^'' 
component, i.e., to a crystallographic C-SHG signal according 
to Table HI Note that out of the six components X^jk^'^ ^^^^ ^re 
allowed x"r^'^ is the only one that was actually observed. Its 
occurrence is limited to the (100) sample where it was sup- 
pressed henceforth by setting the polarization of the incident 
light along y or z, but not in between. 



B. Commensurate SHG contributions 

Figure|3]shows the SHG spectra on a (100) and a (001) sam- 
ple in the multiferroic AF2 phase at 8 K. Rich spectra showing 
all the transitions in Fig. |2a) are observed. We identify Xyyy, 
Xy~, and Xyxx- According to TableUthe observed SHG contri- 
butions reflect the change of point symmetry from 2y/myl' to 
2yl' by the emergence of the spontaneous polarization. The 
assignment is corroborated by temperature dependent mea- 
surements in Fig. |3jc). The ED-SHG contributions in Fig. [3] 
are non-zero only when P^^ ^ 0. All contributions behave ac- 
cording to xfj'i( ^ 

which, because of P^ oc (r^^, _ r)i/2 

and /sHG °^ l^v^P' leads to a linear temperature dependence of 
the SHG intensity in the AF2 phase. 

We thus see that in spite of the IC magnetic order of 
MnW04 the SHG signal can be fully understood on the ba- 
sis of the macroscopic point group symmetries thus far. The 
SHG signal picks up the magnetically induced spontaneous 



net polarization as commensurate projection of the IC spin or- 
der while SHG contributions coupling to the underlying IC 
magnetic structure are not detected. Neither the polarization- 
dependent nor the temperature-dependent SHG measurements 
exhibit peculiarities revealing the presence of an aperiodic 
phase. 

Another noteworthy issue in Fig.[3]!b) is the exceptionally 
pronounced resonance of the component which is due to 
phase matching4i Nearly noncritical type I phase matching 
can be realized in MnW04 for incidence along x at 1 .95 eV, 
as can be seen from the refractive indices of the crystal.— 
This constitutes a rare example of phase matching activated 
by magnetic order and, because of the associated spontaneous 
polarization, presumably the first occurrence ever of "multi- 
ferroic phase matching". 



C. Incommensurate SHG contributions 

Figure|4]shows the SHG spectrum obtained on a (010) sam- 
ple in the multiferroic AF2 phase at 8 K with k\\y. A SHG 
signal from the ^Aig "^Tig transition with a variety of most 
unusual properties is observed. First, SHG of the ED, MD, 
and EQ type is forbidden according to Table U if the macro- 
scopic point symmetries 2y/myl' or 2,1' are applied. Second, 
the temperature dependent SHG data reveal that the SHG sig- 
nal emerges in the non-multiferroic AF3 phase (with a change 
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of slope at Tap:) in contrast to the SHG signals in Fig. Oc). 
Third, the polarization dependent data in Figs. Htb) andHfc) 
reveals lobes consistent with a coordinate system that is ro- 
tated by about 25° around the y axis of the Cartesian system. 
Because of the violation of the inversion symmetry in the AF2 
phase, the observed SHG contribution can be assigned to an 
ED process. The corresponding susceptibility is X^J^u if the 
rotated coordinate system (.x', y' , z') shown in Fig. |4] is ap- 
plied. Note that the axes x' and z' do not only deviate from the 
X and y axes but also from magnetic easy or hard axisi^ We 
speculate that the rotation may be related to the propagation 
direction perpendicular to the planes defined by the IC prop- 
agation vector k^F,. This direction includes an angle of 28° 
with the X axis. All the aforementioned inconsistencies sug- 
gest to associate the "forbidden" SHG signal in Fig.|4]to the 
IC magnetic order which has not yet been included into the 
symmetry analysis. This will be discussed in the following. 

Remember that for a symmetry analysis as in Table 1 only 
global symmetries are considered. Site symmetries, defects, 
and translations, i.e., local properties that are spatially con- 
fined to the extension of the unit cell, are neglected so that 
only the point group symmetry of a crystal is considered. For 
SHG experiments in the optical range this is a reasonable ap- 
proximation because the wavelength of the probe light ex- 
ceeds the extension of the unit cell by three orders of magni- 
tude. However, IC phases are characterized by structural mod- 
ulations exceeding the extension of the unit cell. Therefore the 
spatial variation of the electromagnetic light field across one 
period of the IC wave vector may no longer be negligible so 
that the light begins to sense the local symmetry of the IC 
structure!^ 

Based on the Neumann principle observation of SHG from 
X'fjj indicates the presence of a local symmetry lower than 
2vl'. Here, the only matching point group is 11'. This be- 
comes reasonable once the local structure of MnW04 is con- 
sidered. The spatial variation of the electromagnetic light field 
couples to the local violation of the twofold rotation symmetry 
2v due to the IC nature of the ellipsoidal magnetic spin wave 
and leads to the reduction of the local symmetry from 2,1' to 
11'. 

Figure |5] illustrates such a violation of the local symmetry 
by IC order. A section whose extension is determined by the 
coherence length of the SHG process is considered so that any 
SHG light generated in the crystal section contributes coher- 
ently to the SHG yield. Figure 13 a) shows a commensurably 
modulated chain of atoms. One can see that the commensu- 
rate modulation conserves the symmetry axis 2,, that is indi- 
cated by the twofold rotation axis in the sketch. In contrast, 
the IC modulation in Figure |3b) violates the lattice periodic- 
ity and destroys the rotational symmetry within the considered 
section. Hence, locally, i.e., in the finite section probed coher- 
ently by the SHG process, the 2,, symmetry is lost. 



V. CONCLUSION 

In conclusion, we have analyzed the coupling of optical 
SHG to an incommensurately ordered state using multiferroic 
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FIG. 5: Model explaining the presence of "forbidden" SHG contri- 
butions in crystals with IC order, (a) A commensurate modulation 
of the crystal structure retains the twofold axis of the structure, (b) 
An IC modulation violates the twofold symmetry locally. Additional 
SHG contributions may therefore emerge. 



MnW04 as model compound. Two principally different SHG 
contributions coupling to the IC order parameters were iden- 
tified. The first one couples to a commensurate projection of 
the IC state and is typically represented by a secondary, in- 
duced order parameter. This SHG contribution can be under- 
stood on the basis of a point symmetry analysis for commen- 
surate structures and is indistinguishable from SHG contri- 
butions induced by primary commensurate order obeying the 
same point symmetry. The second SHG contribution couples 
to the primary IC order. It reflects the local symmetry of the 
crystal on the length scale of the IC modulation. Since the 
IC order breaks local symmetries, additional SHG contribu- 
tions that are forbidden by the global point group symmetry 
of the infinitely extended crystal can emerge. They uniquely 
identify the IC state and can be used to image IC domains and 
separate them from any coexisting commensurate domains. 
Here, IC domains expand the established concept of ferroic 
domains because in contrast to commensurate translation do- 
mains neighboring IC domains differ by a lattice translation 
that exceeds the boundary of a single unit celli^ 

From the point of view of application we have character- 
ized SHG as a tool for studying the properties of coexisting 
commensurate and IC forms of order in a single experiment. 
The local symmetry reduction in an IC crystal leads directly 
to the manifestation of additional, "hidden" degrees of free- 
dom for forming domains that are accessible by SHG. These 
additional degrees of freedom play an important role for the 
manifestation of the functionalities of materials with IC order. 
For instance, in multiferroic ]VlnW04 the correlation between 
the magnetically ordered IC state and its domains on the one 
hand and the commensurate dielectric properties of the mate- 
rial on the other hand is a key for understanding the complex 
magnetoelectric interactions of the crystal. 
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